Model of the electron-phonon interaction and optical conductivity of Bai_ x Ka;BiC>3 superconductors 
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We investigate the physical properties of the Bai-^K^BiOa compounds with a focus on the optical proper- 
ties. Results from the simple Holstein model, describing a single band coupled to an oxygen breathing mode 
with parameters derived from first principles calculations, are in excellent agreement with a broad range of ex- 
perimental information. It accounts for an insulating parent compound at x — with a direct- (optical) and 
an indirect-gap, and a metal insulator transition around x = 0.38. Strong electron-phonon coupling leads to 
spectral weight redistribution over a frequency scale much larger than the characteristic phonon frequency and 
to strongly anharmonic phonons. We find that the metallic phase in the vicinity of phase boundary is close to 
the polaronic regime. 

PACS numbers: 74.25.Gz, 74.70.-b, 71.38.-k 



Compounds of the Ba1_3.K3.BiO3 family exhibit both su- 
perconductivity (T c w 30 K) and an exotic semiconducting 
state [1-5]. The parent compound (BaBiOs), with one va- 
lence electron per Bi atom, is an insulator. Upon substitution 
of Ba with K an insulator-metal transition occurs at x « 0.37. 
Significant effort has been made to understand optical con- 
ductivity (OC) measurements, which, in the metallic phase, 
display in addition to a Drude peak a substantial amount of 
spectral weight in the mid-infrared (MIR) frequency region 
[6, 7]. However, there is no consensus on the origin of the 
features of the OC. Interpretations have been put forward ar- 
guing in favor of very weak e-ph coupling [8-10] as well as 
very strong e-ph coupling, which can even be responsible for 
polaronic and bipolaronic behavior [5, 11-15]. 

In this letter we reevaluate the physical consequences of the 
e-ph interaction in the Ba1_3.K3.BiO3 family, within a simple 
model with parameters extracted from recent electronic struc- 
ture calculations [16, 17]. To solve the model we use Dynam- 
ical Mean Field Theory (DMFT) [18], with an Exact Diago- 
nalization solver. We focus on the impact of e-ph coupling 
on the OC, and use comparison of theory and experiments to 
elucidate the strength of the e-ph coupling in this class of ma- 
terials. 

Band theory calculations indicate that a single anti-bonding 
Bi-O sp(a) band crosses the Fermi level [19, 20]. Further- 
more, only optical phonons couple strongly to the electrons. 
Our central approximation is to retain the coupling to a sin- 
gle phonon mode, the breathing mode, and take its coupling 
constant to be independent of momentum (this derivation is 
described in the supplemental materials (SM)). This mode is 
most important since its condensation gives rise to the charge 
ordered insulating state at low filling [21]. Therefore we use 
a minimal single band model Hamiltonian with Holstein e-ph 
coupling (see SM). 

A quantity of central importance is the effective e-ph cou- 
pling strength defined by the low frequency behavior of the 
electron self-energy. In DMFT, the effective mass of the 
electron is given by nih/m* = (1 — 9ReS(w)/9w|„ = o) , 
which is related to the effective coupling, A e ff, by m*/rrib = 
1 + A c ff. mb is the electronic band mass of the carriers. The 



value of A c ff includes phonon renormalization effects, and will 
generally be significantly higher than the bare coupling, Ao, 
defined as (2g 2 /hu!o)N(eF), where g is e-ph coupling con- 
stant, wo is phonon frequency and N(ep) is the density of 
states at the Fermi level. A second way to estimate the ef- 
fective or renormalized value of lambda is via the formula, 
A off = ~g 2 N{e F )RcD{0) = A /(l + 2R C n(0)/fto;o), where 
D(uS) denotes phonon propagator and n(w) is phonon self- 
energy. In our case the above definitions lead to the same 
value for A e ff. It is the renormalized value of the e-ph cou- 
pling, i.e. A c ff, that should be compared with experiments 
[22, 23]. 

In DMFT [18], due to the absence of vertex corrections, the 
OC can be expressed as a functional of the fully interacting 
single particle Green function, G(_k, hoj), as 
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where Ak a (e) = ( — l/7r)ImG cr (ek, e) is the interacting spec- 
tral function, v a (k) = ^p-, /(e) is the Fermi function, and 
Vcell = Nbio, 3 is the unit cell volume. Nbi is the number of Bi 
ions in a unit cell. The lattice constant value is a = 4.27 A. A 
Landau parameter F\ [24] renormalizes the current operator 
inEq. (1) [25]. 

When only intra-band optical transitions relative to the low- 
est conduction band contribute to the optical spectral weight 
(OSW), one obtains the restricted or partial sum rule [26] 
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From a knowledge of the optical spectral weight, we can cal- 
culate the plasma frequency, given by W op t = u 2 /4. 

To appreciate the range of frequencies over which the opti- 
cal spectral weight redistributes, it is customary to define the 
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effective carrier number per Bi ions participating in optical 
transitions, N c g(uj), defined by the partial integral of the OC 
for frequencies less than cj as 
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where ?no is the free-electron mass. In order to compare the 
theoretical results with the experiments of Ref. [7] for N e g (ui) 
we adopt their convention and set too / m b = 1. 

Integrating up to infinite frequency one recovers the OSW. 
A reduction of W op t relative to its non-interacting value, as 
given, for example, in a LDA estimation [27], is usually taken 
as a sign of the presence of electron-electron correlation [28]. 
Before focusing on the Bai_ x K x Bi03 compounds, we now 
revisit this issue on general grounds for the e-ph problem. 

Fig. 1 (a) displays the normalized OSW in the normal and 
CDW phases at half-filling (x = 0) for a particle-hole sym- 
metric model. The most important feature of Fig. 1 (a) is that 
in the weak to intermediate range of Ao the total OSW in a ho- 
mogeneous normal metallic phase depends only weakly on the 
strength of the e-ph coupling as in Migdal-Eliashberg theory. 
The OSW reduction due to e-ph coupling is only appreciable 
by approaching the polaronic regime, but even upon entering 
that regime the OSW is substantial while the effective mass is 
considerably enhanced, as shown in Fig. 1 (a) and (b). These 
results are similar away from half-filling. Fig. 1 (b) shows the 
Drude part [29] normalized by the non-interacting total optical 
integral, W op t(g = 0). In our case, without any other interac- 
tion, W opt (g = 0) = W D (g = 0). The W D {g)/W D {g = 0) 
and inverse effective mass agree and both go continuously to 
zero at the critical coupling strength, Ao c ~ 0.455. 

Hence the variation in the Drude weight or the effective 
mass with e-ph coupling is more prominent than the variation 
in OSW. Fig. 1 (a) also shows that for Ao < Ao c , ordering 
of itinerant carriers in a charge density wave (CDW) phase 
leads to a larger reduction of OSW than normal phase. How- 
ever, for Ao > Aoc, the carriers are localized in real space, due 
to bipolaron formation, and CDW ordering enhances OSW. 
Therefore, the e-ph coupling dependence of the difference be- 
tween OSW in the metallic and CDW phases is not monotonic 
(similar to the temperature dependence of the OSW [31]). 

To focus on the Bai-^K^BiOs compounds, we deter- 
mine the hopping parameters to fit the LDA band structure 
[19] which results in following hopping parameters t\ = 
0.3926, t 2 = 0.0516, t 3 = -0.0017, U = -0.0987, all in 
eV. The strong energy variation of the electron DOS for the 
real material (see SM) leads to a strong doping dependence of 
the bare e-ph coupling, Ao, which plays a central role for un- 
derstanding the physical properties of the Bai_ K K x Bi03 fam- 
ily. For the phonon frequency and e-ph coupling we choose a 
set of the bare parameters, Hluq = 80 meV, and g = 0.2275 
eV. The value of the e-ph coupling to the breathing mode, g, 
matches closely to the recent results of first principles calcu- 
lation of the deformation potential at the zone boundary [17] 
(see SM). The characteristic phonon frequency represents a 
bare parameter. Below, we shall see that it becomes renor- 
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FIG. 1. (Color online) Panel(a): normalized interacting OSW as 
a function of bare e-ph coupling, Ao, at half-filling (x = 0) for 
a particle-hole symmetric model in the normal and CDW phases. 
Panel(b): normalized Drude weight and inverse effective mass as a 
function of Ao. The TB parameters for this figure are: t\ — 0.3926 
eV and &j = 0.08 eV. 



malized, so that the final value is representative of the exper- 
imental value. For simplicity we take all the bare parameters 
in the Hamiltonian to be independent of doping (rigid band 
picture). Considering the simplifications made, our results are 
in surprisingly good agreement with experiments. 

For the parent compound, BaBiOa (x = 0), we show the 
calculated conductivity as a function of frequency in panel 
(a) of Fig. 2, along with two measurements [6, 7] performed 
at room temperature. The calculated conductivity displays a 
characteristic insulating spectrum; the Drude peak is absent, 
and spectral weight has been transferred to the region above 
the gap. The calculated gap value, A opt ~ 2 eV, and the 
subsequent absorption strength are similar in magnitude to the 
observed values [6, 7]. The inset shows the phonon density 
of states in CDW phase with a renormalized peak near w 72 
meV, in good agreement with the experimental value of the 
breathing mode phonon frequency (Experimental data in [20]) 
[33]. 

In Fig. 2 (b) we show N e g(uj) vs. uj for the undoped sys- 
tem. The calculated iV c ff is in good agreement with the exper- 
imental data in the range of frequencies in which inter-band 
transitions play no role (u> < 3 eV). Fig. 2 (c) shows the elec- 
tron DOS in the charge ordered phase; it shows a gap value 
ss 0.5eV, which corresponds to the indirect gap in BaBiOs 
[34]. 

Upon doping the system we find that at doping value around 
x w 0.38, charges distribute uniformly between the two sub- 
lattices. Investigation of the electronic DOS indicates that this 
is a transition to a metallic state with a finite DOS at the Fermi 
energy. In the intermediate doping range < x < 0.38, the 
electron density oscillates between two values in successive 
iterations of the DMFT cycle and a converged answer was 
not found. This implies that the system either has a tendency 
towards another phase, like a non-commensurate CDW, or a 
tendency towards phase separation, both of which cannot be 
seen in our study [35]. 

Fig. 3 (a) shows the calculated OC for two doping val- 
ues in the metallic phase. As the system is doped into the 
metallic phase, the optical response shows a Drude peak, cen- 
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FIG. 2. (Color online) (a) Calculated and measured OC at half-filling 
(x — 0). Measured OCs of BaBiC>3 are reproduced from [6, 7], la- 
beled by (1) and (2), respectively. Inset: The renormalized phonon 
density of states, (b) Calculated and measured effective carrier num- 
ber per Bi ion participating in the optical transitions, (c) Electron 
DOS in the charge ordered phase with an indirect gap of the order of 
A ~ 0.5 eV. A and B indicate the majority and minority sublattices 
for electron density. 



tered at ui = 0. This indicates coherent quasi-particle motion, 
and an additional mid-infrared (MIR) feature arises, caused 
by single- and multi-phonon excitations and absorption pro- 
cesses. The MIR contribution is an incoherent component of 
the OC, and corresponds to optical transitions in which the 
carrier is excited but leaves behind a local distortion, which 
corresponds to real excitations. This MIR component occurs 
at frequencies characteristic of these excitations. As we saw 
in Fig. 1 (b), the e-ph interaction causes a reduction of the free 
carrier Drude peak weight; in fact it is approximately reduced 
by the quasi-particle renormalization factor (1 + A c ff ) _1 and 
the spectral weight increases in the MIR energy range reflect- 
ing the incoherent scattering excitations that result from ab- 
sorption processes assisted by the e-ph interaction. The MIR 
frequency range extends well beyond that of the characteristic 
phonon frequency. 

In the inset of Fig. 3 we show the effective carrier num- 
ber calculated from OC. In the metallic phase N e g(uj) at first 
increases steeply at low frequencies, on account of the appear- 
ance of the Drude band centered at u = 0, and continues to 
increase in the mid-infrared region, eventually saturating at a 
high frequency. The trends are in good agreement with the 
experimental data [7, 36]. 

Fig. 3 (b) shows the phonon density of states. The bare 
phonon frequency is softened and the phonon DOS shows 
structures similar to those observed experimentally [32]. The 
effective oscillator potential and anharmonicity is further ex- 
plored in the SM. Consistent with experimental results [37], 
the effective potential shows a double well structure, with one 
minimum disappearing with increased doping. 

We now turn to the doping dependence of the OSW. Due to 
a higher value of the DOS at (x = 0) the e-ph coupling (Ao) 
is largest in the parent compound and decreases with doping. 
Fig. 3 (c) shows the normalized interacting OSW vs. doping. 
It is clear that the OSW is greatly reduced for small values of 




FIG. 3. (Color online) (a) Calculated and experimental [7] conduc- 
tivity for the indicated doping x at zero temperature. Inset: Effective 
carrier number participating in the OC as a function of frequency, (b) 
The phonon density of states at x = 0.5. (c) Normalized interacting 
OSW by non-interacting value as a function of doping, (d) Drude 
weight, Wo, of the OC, normalized by total weight Wopt, as a func- 
tion of the doping level. In panels (c) and (d), system at x — is in 
CDW phase, while other data points show normal metallic phase. 



x due to the e-ph interaction. In the metallic phase, upon in- 
creasing the doping level, electrons become more undressed, 
which leads to a smaller optical weight reduction. 

Fig. 3 (d) shows the Drude pan (the MIR part is W op t — 
Wd) normalized by the total optical integral as a function of 
doping at T = [38]. In the metallic phase the OSW redis- 
tributes between the Drude peak and the MIR shoulder, and 
by increasing the doping level the Drude weight rises very 
quickly. For instance, the Drude peak weight increased by 
almost a factor of two with increasing x from x = 0.4 to 
x = 0.5. 

In the DMFT framework, the effective e-ph coupling can 
be extracted either from the phonon propagator at zero fre- 
quency, or the effective mass, or from the relative weight of 
the Drude peak to the total OSW in the absence of e-ph cou- 
pling, W B (g)/W opt {g = 0) w (1 + Acff)- 1 . We found that 
the effective e-ph coupling is strongly doping dependent with 
the highest value (A e ff > 4) at x = 0.4 close to the transition 
to the insulating phase. With such a high value of the effective 
e-ph coupling the formation of an inhomogeneous phase in the 
vicinity of the phase boundary is highly likely. At higher dop- 
ing levels where the system is likely to be homogeneous, we 
find A c ff f» 1.6 for x = 0.5. 

The calculated plasma frequency for x = 0.5 is H 2 ujp = 
8.95 eV 2 , again in good agreement with the experimental 
value [15]. Using the calculated values for the plasma fre- 
quency for x = 0.5 we can use the high temperature behavior 
of the measured resistivity to estimate the e-ph coupling, A^ 
[39]. Using the high temperature resistivity formula [14], 

p(T)«p ph (T)«^^T (4) 
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we estimate a value of X y cS ' « 1.75. 

In support of strong e-ph coupling, heat capacity mea- 
surements for Bai_j.Kj.BiO3 w i m x ~ 0.4 show a jump 
in the specific heat of the order AC P (T C )/T C w 4 - 
5 mJ mol _1 K -2 at T c [40]. One can relate this to the co- 
efficient of the electronic specific heat, 7, by the BCS ex- 
pression AC P (T C )/T C = 1.437 [41] ■ Then from the rela- 
tion 7 = (2/?>)it 2 N(e f )fc|(l + Acs) with use of N(e F ) = 
0.23 states/eV spin cell for x = 0.5, A e ff is estimated to be 
1.6 — 2.2, consistent with our previous estimates. 

In conclusion we have revisited the physical properties of 
Bai-jK^BiOs, using a model Hamiltonian with parameters 
inspired by recent first principles calculations [ 16, 17] but go- 
ing beyond these calculations by taking into account the full 
many body physics and inherent anharmonicity of a strongly 
coupled e-ph system. For small doping the coupling is strong 
enough to sustain polaron and small bipolaron formation and 
the coupling decreases for large doping. For the metallic 
regime near the phase transition boundary, we find a situa- 
tion intermediate between a Fermi/BCS liquid and a bipola- 
ronic superconductor, in agreement with experimental stud- 
ies [42, 43]. The model accounts for experimentally deter- 
mined gaps (both direct and indirect) in the parent compound 
at x = and an abrupt insulator transition at a critical dop- 
ing. In the metallic phase, we find a strongly anharmonicity 
in agreement with recent experiments [37] and transfer of sub- 
stantial electron phonon coupling strength to low frequencies 
as shown in the phonon density of states. Strong e-ph cou- 
pling reduces the width of the Drude peak and transfers spec- 
tral weight to the mid-infrared structure. Spectral weight re- 
distributes over a frequency scale that is much larger than the 
phonon frequency. Further extensions will require the incor- 
poration of couplings to other phonon modes and of additional 
bands in the modeling. These are not expected to change dras- 
tically the main conclusions presented in this paper which fo- 
cused on the electronic properties but are certainly required 
for a more realistic study of the phonon spectra. 
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Model of the electron-phonon interaction and 
optical conductivity of Bai_ x K x Bi03 
superconductors - Supplemental materials 



I. ELECTRON-PHONON COUPLING 

Here we provide additional information on the derivation of 
a Holstein-type e-ph coupling from the Rice-Sneddon model 
Hamiltonian [1]. The phonon part and e-ph interaction in the 
Rice-Sneddon model Hamiltonian is given by 



H ph + H e _ ph = ^(«- a + u%u? a ) 

ia 

- Q ni(u ia+ - Ujq-) (SI) 

ia 

where Ui a ± denote displacements of two oxygen ions at 
sites Ri±ae Q /2, placed in the ath direction at a distance a/2, 
and Q is the strength of this interaction. Mq is the oxygen ion 
mass, and u>o is the phonon frequency. 

Introducing the boson representation for each oxygen coor- 
dinate, u ia = y/h/2M o u) (B\ a + B la ), Eq. (SI) reads 



Hph + He-ph — fa^O y^(-B|q-Bm + -) 

ia 

-s{^X nA[Bl+ -^- )+h - c] - (S2) 

ia 

With these phonon coordinates, being attached to oxygen 
ions, there is a nonzero overlap between nearest-neighbor op- 
erators, i.e., they are not orthogonal. It is better to define 
new orthogonal boson operators [2]. By Fourier transform- 
ing the boson operators, Bi a ± = (1/y/N) cxp[ik.(Ri ± 
ae a /2)]B\ (ia ±, it is easy to show 



v i k 

(S3) 

where is the number of sites, 6k represents a symmetric 
phononic mode 



bk = —— ^ BkaSOT^, (S4) 

Mk ^ 2 



and 



H k = Jsin 2 (k x /2) + sin 2 (k v /2) + sin 2 (k z /2). (S5) 



6 



The new operators, b^, fulfill the boson commutation rela- 
tions. By transforming the new operator back to real space 
we obtain the following form of the phonon related terms 



i 

-^SiMtf + bj) (S6) 



where <7y = 2Q h/2Mo^o"fij ■ The coefficients 7^ = 
Sk/^k exp[— ik.(Ri — Rj)] are responsible for the long- 
range e-ph interaction, which is a consequence of orthogo- 
nalization of the local phonon modes. However, these terms 
fall off rapidly with distance (e.g. 7(0) = 1.19, and 7(1) = 
—0.11) [2]. In the Holstein model, we only keep the first term 
g = gu = 2Q ^ h/2MoUo~f(0). Therefore, we use following 
model Hamiltonian: 



+ gJ2(m - l)(bi + b\) + hwv (S7) 

i i 

where ci a (c| CT ) and bi (b\) are the destruction (creation) 
operators for electrons with spin a and local vibrons with fre- 



J 



is the 



quency ujq, respectively. The operator n, = ^ 
electron density at site i, the parameters are the electron 
hopping matrix elements between sites at i and j, and g de- 
notes the e-ph coupling constant. 

Using Hlj^ da = 0.065 eV, and Q = X> LDA /6 the value of 
g would be g = 0.237 eV. Here X> LDA = 13.3 cVA" 1 is the 
deformation potential obtained from the LDA calculation [3]. 

The displacement of the oxygen ion is approximately given 
by u la w (l/3)Vfc/2M wo7(0)(&i + b\). Our model 
Hamiltonian calculation, with Hu!q = 0.08 eV, gives a value 
(bi + b\) w 4.5 for the undoped case, which leads to a breath- 
ing mode distortion of the order of m ~ 0.073 A, in agreement 
with the experimental value 0.085 A [4]. 



H. DENSITY OF STATES 

In Fig. SI we show the resulting density of state (DOS) for 
a specific choice of hopping parameters t\ = 0.3926,^2 = 
0.0516, t 3 = -0.0017, t 4 = -0.0987, all in eV (TBI). These 
correspond to successively further neighbor hopping. Note the 
appearance of a (diverging) van Hove singularity, the result 
of including beyond nearest neighbor hopping, and in strong 
contrast to the rather flat DOS for the nearest neighbor model 
(TB2). The strong energy variation of the electron DOS for 
the real material leads to a strong doping dependence of the 
bare e-ph coupling, Ao, which plays a central role for under- 
standing the physical properties of the Bai-^K^BiOa family. 
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FIG. SI. (Color online) The density of states of the tight-binding 
model (TBI). Note the van Hove singularity. Vertical lines indicate 
the non-interacting chemical potential for x = 0.5, x — 0.4 and 
x = 0. Also shown is the DOS for nearest neighbor hopping only 
(TB2). 





FIG. S2. (Color online) The phonon displacement probability- 
distribution function (left panel) and the effective phonon potential 
(right panel) for three different doping levels. 



m. EFFECTIVE PHONON POTENTIAL 

In an electron-phonon system, an effective phonon poten- 
tial, V e s, can be defined using an effective wavefunction de- 
fined by ip c s(q) = y/P(q), which is taken to be a solution of 
the one dimensional Schrodinger equation [5], 



d 2 ^ cS {q) 
dq 2 



V c s(q)tpcs(q) = £^eff(?)> 



(S8) 



where q = b + b^ and P(q) = (4>o\q) (q\4>o) is the phonon 
displacement probability-distribution function (PDF) [6]. The 
wave function \4>o) is the ground-state wave function and 
\q)(q\ is the projection operator onto the subspace where the 
phonon displacement operator, q, at a given site has value q. 
This quantity is a measure of the distribution of the local dis- 
tortions. Having P(q), the effective phonon potential V c s(q) 
is given by 

VM = E+ —{Tw-2{m) )■ (S9) 

In Fig. S2 we show results for the probability distribu- 
tion function and corresponding effective potential as deduced 
from Eq. (S9) for three different doping levels (we ignore the 
constant E, so we only discuss the variation in the shape 
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of the potential and not its absolute value). In each case 
P{q) has a narrow peak which shifts to slightly larger q with 
increasing doping. The average value of the displacement, 
(q) = J_ dq qP{q), agrees with a direct calculation of 
(6 + 6+) and is given by -(2g 2 /hu )((n) - 1) [5]. The effec- 
tive phonon potential shows an asymmetric double well struc- 
ture, with one of the minima disappearing with increased dop- 
ing. Near the lowest doping value (x = 0.4), the potential ex- 
hibits pronounced anharmonicity, which eventually decreases 
as the doping increases. For the highest doping level, the shal- 
low minimum has disappeared, and fluctuations decrease, in- 
dicated by the presence of a single minimum in the effective 



potential. The enhanced fluctuations near the instability point 
are a signal that an instability is imminent. 
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